Magnetic anisotropy of FePt: effect of lattice distortion and chemical disorder 
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We perform first principles calculations of the magnetocrystalline anisotropy energy in the ve 
Ll FePt samples studied experimentally by Ding et al. [J. App. Phys. 97, 10H303 (2005)]. The 
effect of temperature-induced spin fluctuations is estimated by scaling the MAE down according 
to previous Langevin dynamics simulations. Including chemical disorder as given in experiment, 
the experimental correlation between MAE and lattice mismatch is qualitatively well reproduced. 
Moreover we determine the chemical order parameters that reproduce exactly the experimental 
MAE of each sample. We conclude that the MAE is determined by the chemical disorder rather 
than by lattice distortion. 

PACS numbers: 75.30.Gw 75.50.Ss 71.15.Mb 71.15.Rf 



Due to its extraordinarily high magnetocrystalline 
anisotropy energy (MAE), Llo FePt is of considerable in- 
terest to the development of ultrahigh density magnetic 
recording applications and spintronics devices. From a 
theoretical point of view, there is an obvious need for 
a complete first principles model of FePt to be used in 
generating effective spin Hamiltonians for the purpose 
of atomistic and multiscale modelling. Amongst many 
other issues, this requires an understanding of the role of 
interfacial effects and chemical disorder. The large effect 
of chemical disorder on the MAE of FePt has already 
been outlined both experimentally^ and theoretically^. 

Recently, the experiments were extended to thin films 
of FePt deposited on different substrates^. A strong cor- 
relation was revealed between the MAE of the FePt sam- 
ple and the lattice mismatch of the FePt films with re- 
spect to the substrate^. The experimental data are sum- 
marised in Table [IJ The chemical order parameter, s, is 
defined as the probability of finding an Fe atom on a nom- 
inal Fe site or, equivalently, as the probability of finding 
a Pt atom on a nominal Pt site. In the experiment, the 
chemical order parameters were derived from the X-ray 
diffraction intensities 7(001) and 7(002), ((xyz) denot- 
ing the plane of diffraction), through the relationship 1 4 
s ~ \J 7(001)/7(002) and normalizing s to unity for sam- 
ple no. 3. We refer to the experimentally obtained chemi- 
cal order parameters as s e for distinction from the chemi- 
cal order parameters s obtained later by fitting calculated 
MAE- values to experiment. 

The aim of the present work is to investigate in de- 
tail the effect of lattice distortion and chemical order on 
the MAE of FePt by means of the relativistic Korringa- 
Kohn-RostokerS"^ method as combined with the coherent 
potential approximatiorPSl (KKR-CPA). In order to dif- 
ferentiate between the two main properties characterizing 
the samples, namely, the lattice distortion and the chem- 
ical disorder, we perform calculations with and without 
the inclusion of chemical disorder. We then fit the calcu- 



TABLE I: Summary of experimental results by Ding et al.^ 
lattice parameters a and c, chemical order parameter, s, mag- 
netocrystalline anisotropy energy per formula unit, K, and 
diffraction intensity ratio, 7(001)// (002). 



Sample 


a (A) 


c(A) 


Se 


K (meV) 


7(001)/7(002) 


1 


3.88673 


3.69977 


0.709 


0.493 


1 


2 


3.88279 


3.69387 


0.978 


0.696 


1.9 


3 


3.89752 


3.68964 


1.000 


0.841 


1.985 


4 


3.89646 


3.69175 


0.965 


0.788 


1.85 


5 


3.86954 


3.71378 


0.615 


0.271 


0.7536 



lated MAE to the experimental values using the chemical 
order parameter, s, as a fitting parameter and draw con- 
clusions from the results of our calculations. We find that 
chemical disorder of each sample is the more important 
factor in determining the experimental 3 MAE. 

As the relativistic KKR method is well documented in 
the literature (see e.g. [7]), here we merely describe some 
details of our calculations. We used Density Functional 
Theory within the Local Spin-Density Approximation 
(LSDA) as parametrised by Vosko et o/W The effec- 
tive potentials and fields were treated within the atomic 
sphere approximation (ASA). As the thin-film samples in 
the experiment had a thickness of approximately 20 nm 
(60 formula units) P surface contributions to the MAE 
should be negligibly small. We therefore modelled the 
FePt samples as face-centered-tetragonal (fct) bulk lat- 
tices with lattice constants as displayed in Table [I] The 
self-consistent calculations were performed by using the 
scalar-relativistic approximation, i.e., by neglecting spin- 
orbit coupling^ and solving the Kohn-Sham-Dirac equa- 
tion using a spherical wave expansion up to an angular 
momentum quantum number of t = 3. As in earlier the- 
oretical workpl we used the coherent potential approx- 
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imation (CPA) to elucidate long-range chemical disor- 
der effects in FePt. In combination with KKR, the CPA 
has proved particularly useful in calculating the physical 
properties of chemically disordered alloys^. The partially 
disordered FePt alloy is modelled by a stack of alternat- 
ing layers with the chemical compositions of Fe s Pti_ s 
and Pt s Fei_ s . 

The MAE is then evaluated using the magnetic force 
theorem^, which states that the difference in a system's 
total energy for two different directions of magnetiza- 
tion can be approximated by the corresponding difference 
of the band energies, neglecting further self-consistency 
i.e., keeping the effective potentials and fields fixed. From 
previous experience we know that for transition metal 
systems these potentials and fields can safely be taken 
from self-consistent scalar-relativistic calculations^. In 
order to achieve a relative accuracy within 5 % for the 
MAE, the associated energy integration was performed 
by sampling 20 energy points along a semi-circular con- 
tour in the upper complex half-plane. At the energy point 
closest to the real axis the ^-integration was calculated 
using 5050 /c-points in the irreducible segment of the two- 
dimensional Brillouin zone. 

As the MAE should vanish at the Curie temperature, 
it is a rapidly decreasing function of temperature. Whilst 
the temperature dependence of the MAE of ordered FePt 
has been previo usly c alculated in terms of different the- 
oretical method j 13 * 14 !, in the present work we do not 
make an attempt to carry out a similar process, since 
site-resolved information is currently not available for a 
chemically disordered system. Instead, for an approx- 
imate comparison with experiments at room tempera- 
ture, we use the scaling obtained for perfectly ordered 
Llo FePt in terms of Langevin dynamics simulations^, 
namely, K T=2 g 3 K ~ 0-6K t =ok- 

Using the methods described above we performed sys- 
tematic calculations of the MAE of each of the FePt sam- 
ples in Table [T] In order to separate the effects of the 
lattice distortion and the chemical disorder, we split our 
study into three stages. In our first set of calculations, the 
FePt samples were modelled as perfectly ordered alloys 
with lattice parameters according to Table [T] As can be 
inferred from Fig. [T] our calculated values spread around 
3 meV/f.u. and show a very minor dependence on the 
variation of the lattice parameters. Moreover, this mod- 
erate variation between the samples is contrary to the 
experimentally observed trend. 

Although high in comparison to experiment, our cal- 
culated MAE values are in good agreement with other 
theoretical results based on the LSDA or the LSDA+U 
approach^. One obvious reason for the discrepancy 
between the theoretical and experimental values is the 
strong temperature dependence of the MAE. We estimate 
this contribution by scaling the calculated MAE down 
by an approximate factor of 0.6, as described above. The 
corresponding MAE- values (also shown in Fig. [I]) are still 
too high as compared to experiment. Thus we conclude 
that, even when taking temperature-induced spin fluctu- 
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Sample No 

FIG. 1: Crosses (solid line): Calculated MAE per formula unit 
for each of the FePt samples in Table [I] modelled as perfectly 
ordered alloys. Circles (dashed line): The same values scaled 
down by a factor of 0.6 in order to account for temperature 
induced effects. Stars (dotted line): The experimental values. 



ations into account, lattice distortion alone can explain 
neither the size nor the trend of the MAE obtained in 
the experiment. 

Subsequently, the chemical disorder of each sample as 
given in Table [I] was taken into account using the co- 
herent potential approximation. The corresponding re- 
sults are shown in Fig. [2] In accordance with earlier 
worlP, long-range chemical disorder drastically reduced 
the MAE; for sample no. 1 (s e = 0.709) we obtained a 
value of 0.4 meV/f.u., while for sample no. 5 (s e = 0.615) 
the MAE almost vanished. In fact, reducing s to 0.5 can 
even cause a change of sign of the MAE. In contrast, for 
samples no. 2 and 4 with a high degree of chemical order 
the MAE was reduced by less than 10 %, and for sample 
no. 3 (s e — 1) the MAE remained unchanged with respect 
to our previous calculations. Taking into account again a 
reduction by a factor of 0.6 due to temperature effects, it 
is obvious that the inclusion of chemical disorder has sig- 
nificantly improved the agreement between experiment 
and theory: the trend of the MAE between the different 
samples is now correct and the magnitudes of the MAE 
are closer to the range reported by the experiment. 
As mentioned above, the chemical order parameters 
we re derived from measured diffraction inten- 
4^ However, due to an incomplete rocking 



in Table I 
sity ratio; 

curved, the measured diffraction intensities, and thereby 
the experimentally obtained chemical disorder parame- 
ters, can only be considered approximate values. Fur- 
thermore, we note the assumption that the sample with 
highest MAE, sample no. 3, refers to perfect chemical 
order, s e = 1. This seems a reasonable working hypoth- 
esis, but one worth investigating theoretically since it is 
central to the interpretation. 

The above uncertainties motivated us to perform a 
third set of calculations, in which the theoretical MAE 
was fitted to the experimental MAE using the chemical 
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FIG. 2: Crosses (solid line): Calculated MAE per formula unit 
for each of the FePt samples in Table [I] modelled as partially 
disordered alloys with the degree of disorder given by the 
experiment. Circles (dashed line): The same values scaled 
down by a factor of 0.6 in order to account for temperature 
induced effects. The experimental values are also displayed 
by stars (dotted line). 
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FIG. 3: Magnetic anisotropy energy (MAE) calculated as a 
function of chemical order parameter, s, for the FePt samples: 
1 +, 2 *, 3 ■. 4 x, 5 •. Solid lines serve as a guide for the 
eye. Open circles are placed at the best-fit chemical order 
parameter for each of the samples. Dashed line: Linear fit. 

order parameter, s, as a fitting parameter. In Fig. [3j for 
each of the samples we present the calculated MAE for 



an appropriate set of chemical order parameters. Firstly, 
for a given sample, i.e. for fixed lattice parameters, the 
theoretical MAE shows a non-linear dependence on s. In 
Fig. [3] the circles indicate the intersection of the calcu- 
lations with the experimental values for each sample as 
indicated. This determines the best-fit order parameter 
that corresponds to the experimental MAE value. As 
can be clearly inferred from Fig. [2j for samples no. 2, 3 
and 4, a smaller degree of chemical order was fitted than 
predicted by the experiment, namely, s ~ 0.836, 0.874 
and 0.863, respectively. In contrast, for samples no. 1 
and 5 an increased degree of chemical order, s ~ 0.782 
and 0.720, was obtained. Although for a given sample 
the theoretical MAE shows a non-linear dependence on 
s, there is a nearly perfect linear correlation between the 
experimental MAE and the best-fit chemical order pa- 
rameters as indicated by the dashed line in Fig. [3j Ob- 
viously, this remarkable linear behavior is the result of 
a subtle interplay of the dependence of the MAE on the 
lattice distortion and the chemical disorder. This is prob- 
ably specific to the data set investigated here rather than 
being a general property. 

In summary, our first principles calculations imply that 
lattice distortion in the FePt samples has only a minor 
effect on the MAE, even opposite to the experimental 
trend. Calculating the MAE using the highly approx- 
imate experimental chemical order parameters signifi- 
cantly improves the agreement between theory and ex- 
periment, in particular with regards to the relative dif- 
ferences in the MAE between the samples. This indicates 
that the substrate-sample lattice mismatch effect on the 
MAE reported by Ding et alP is mainly due to the varia- 
tion in chemical disorder. To circumvent the uncertainty 
of the experimental determination of chemical disorder, 
we, furthermore, determined theoretical chemical order 
parameters that reproduced the experimental MAE val- 
ues. Interestingly, a linear correlation between the MAE 
and the best-fit chemical order parameters is found. It 
should be mentioned that work is underway to perform 
constrained Monte-Carlo simulations of K(T) for chem- 
ically disordered FePt, since this is clearly an important 
factor in relation to experimental data. 
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